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We demonstrate the presence of an important anisotropic magnetoresistance contribution to the 
domain wall resistance recently measured in thin-film (Ga,Mn)As with in-plane magnetic anisotropy. 
Analytic results for simple domain wall orientations supplemented by numerical results for more 
general cases show this previously omitted contribution can largely explain the observed negative 
resistance. 
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The electrical resistance associated with current flow 
across a domain wall (DW) separating uniformly magne- 
tised regions in a magnetic material has been the subject 
of investigation since the 1930s. A large volume of of- 
ten contradictory results exists, mostly attributable to 
the difficulty in separating normally small intrinsic ef- 
fects from the myriad of extrinsic effects that also con- 
tribute in magnetoresistance measurements In recent 
years some consensus has emerged in the study of DWs 
in metallic epitaxial films and nanostructures? where it 
is possible to more fully characterise and control the 
magnetic microstructure, with many results consistent 
with the spin-mistracking models of DW resistance;^ 
Although the small magnitude of the DW resistance lim- 
its potential applications, significant enhancements at 
nanoconstrictions^ allied with advances in the atomic- 
scale control of materials, raise hopes of practical mag- 
netoresistive devices, whilst a close relationship with 
the phenomenon of current-induced DW motion^i^ also 
makes the understanding of DW resistance of consider- 
able importance. 

It has been noted that enhanced magnetoresistance 
effects may be associated with DWs in ferromagnetic 
semiconductors due to the longer Fermi wavelength^ and 
large exchange splitting relative to band widthi^ Riister 
et al^ have observed an 8% increase in the magnetore- 
sistance due to DWs pinned at < 10 nni constrictions 
in in-plane magnetised (Ga,Mn)As nanostructures, and 
Chiba et al}*^ report a significant positive DW magne- 
toresistance in perpendicularly magnetised (Ga,Mn)As 
layers consistent with the theory of Levy and Zhangii 
Tang et al<^ have studied the resistance of 30 — 100 ^m 
devices patterned from in-plane magnetised (Ga,Mn)As 
epilayers. By measuring the average resistance (i?) along 
the sides of the device channel as a 90° domain wall is 
driven through by current pulses, they find a small resis- 
tance drop. Scaling by a wall width of 10 nm. Tang et al. 
deduce a DW resistivity as large as Agi/p —100% — a 
remarkable result implying resistance free current trans- 
port through the region occupied by the DW. Although 
theories exist that predict a negative intrinsic DW resis- 
tance, either as a result of modifications to quantum in- 
terference phenomena^'^ or differences in spin-dependent 
relaxation rateSj^'^ this result is many orders of magni- 



tude greater than any negative DW resistivity previously 
reported in a metal. ^ 

In this work we demonstrate the existence of a size- 
able extrinsic contribution to the negative DW resistance 
measured using the experimental configuration employed 
in Ref. llj. This anisotropic magnetoresistance (AMR) 
effect arises from the circulating currents induced by the 
abrupt change in the off-diagonal resistivity at 90° DW, 
and persists even after the resistance is averaged across 
the sample. An analytic expression is derived for the 
simplest DW orientation, supplemented by numerical re- 
sults for the more general case, which also allow us to 
simulate the experiments where a current-driven DW is 
moved through a microdevice. 

We consider the current flow within an infinitely long 
thin conducting sample with rectangular cross section, 
width w, thickness t. The sample lies parallel to the xy 
plane with the long edge parallel to the x axis (Fig. [1]). 
A dc electrical current / flows through the sample. Ideal 
probes are attached and measure the potential at points 
on either side of the sample separated by a distance I. 

The current density J{x, y) is assumed to be uniform 
as X ±cx), and no current flows through the sides of 
the device: 

J(±(X3,2/) = (/M, 0) = (.7, 0) (la) 
Jy{x,{)) = Jy{x,w) = 0. (lb) 

Within the sample, J is found by satisfying current con- 
tinuity and the steady state Maxwell equation, with the 
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FIG. 1: The geometry considered in this work. A domain 
wall inclined at an angle i? separates regions in which the 
in-plane magnetisation lies at angles ipi and (p2- Potential 
measurements are made using probes A, B, C and D. 
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electric field E and current density related via Ohm's law 



V-J = 

V X E = 
E = pj 



(2a) 
(2b) 
(2c) 



with p a spatially varying resistivity tensor. With in- 
plane magnetised (Ga,Mn)As epilayers, the resistivity in 
directions parallel (py) and perpendicular {pi_) to the 
magnetisation differ^^ with p|| < pj.. If </? is the mag- 
netisation direction in a given domain (Fig. [T]), the cor- 
responding cartesian tensor is 
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^ cos 2lp ^ sin 2lp 
^ ' f sin 2(p 1 - I cos 2ip 



(3) 



where p = {p±_ + p\\)/2, (3 = {p\\ - pi_)/p. With cu- 
bic anisotropy, domain walls divide regions which have 
magnetisation directions differing by 90°, but in the ex- 
periments recently reported this is modified by a weak 
in-plane uniaxial anisotropy. We denote by the angle 
which the normal to the wall (mdw) makes with respect 
to the X axis. Fig. [1] shows the geometry in the case of a 
single domain wall within the device. 

At the wall itself, the boundary conditions are conti- 
nuity in the normal component of the current and in the 
tangential component of the electric field: 



J L nDw 
E II now 



Jx cos d + Jy sin (4a) 
—Ex sin -d + Ey cos -d 

-{QxxJx + QxyJy) siu I? 

+ {QyxJx + QyyJy) COSd. (4b) 



Matching (|4bp when the resisitivity tensor elements are 
different on either side of the DW is not possible with 
a uniform current J = (j,0); the DW induces circulat- 
ing currents and it is these that give rise to an AMR 
contribution to the resistance across the wall. 

To see this, we first consider the case where the de- 
vice channel contains a DW in the yz plane — 0) at 
X = xq, separating regions in which the in-plane mag- 
netisation is at an angle ipi in region 1 and <p2 in region 
2. The average of the longitudinal resistances measured 
along opposite sides of the device, used in Ref. to 



eliminate contributions from the planar Hall effect, can 
be expressed in terms of a difference in voltages at either 
end of the device channel: 

-I{R) = iVB-VA)/2 + iVD-Vc)/2 

- {Vb + Vd)/2-{Va + Vc)/2. (5) 

(The minus sign is because the potential falls along the 
direction of positive current flow.) We take the length 
I (and xo) to be large enough so that the static eddy 
currents induced by the DW are fully contained within 



the area of the device defined by the 4 probes, and then 
the current has its asymptotic value at both x = and 
X — I. This means that the electric field in the y-direction 
is constant both between C and A, ^^^(O, y) = Py^h ^nd 
D and B, Ey{l,y) = Py^j, and the voltage changes lin- 
early between points C and A, and between D and B. 
The voltage averages in ([5]) can then be re-expressed as 
integrals of the voltage across the device, e.g. 



{Vb+Vd)/2^- [ V{l,y)dy, 
w Jo 



(6) 



and the two terms in ([5]) then can be combined using 

V{1, y) — F(0, y) ~ — /q Ex{x, y)dx to give (i?) in terms 
of an integral of the electric field over the area of the 
device between the probes: 



I{R) = - [ [ Ex{x,y)dxdy. 

w Jx=a Jy=0 



(7) 



Splitting the integral into separate contributions from the 
two domains within each of which the resistivity tensor 
is constant, and using the following results 



Jx{x,y)dy ^ jw, 



Jy(x, y)dx = (8) 



that are found by integrating the continuity equation over 
regions Vl = {{x',y) e R'^ \ < x' < x,0 < y < w} and 



n = {ix,y') 
yields 



< X < 1,0 < y' < y} respectively. 
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x=0 Jy=0 



Jy{x,y)dxdy. (9) 



The first term describes a resistance that linearly interpo- 
lates between the asymptotic resistances of the channel 
in the two uniform magnetisation states. In Ref. (llj 
differences between measured resistance values and this 
linear interpolation have been interpreted as originating 
from an intrinsic DW resistivity. However, the final term 
in ([9]), which henceforth we denote i?AMR, is a new con- 
tribution that we find, which results directly from the 
discontinuity in the resisitivity at the DW and which is 
proportional to the total parallel current induced on ei- 
ther side of the DW. 

To obtain an estimate for the value of this additional 
contribution to the DW resistance, we consider the case 
where Lpi = —ip2 — f which applies when the hard axis 
is perfectly aligned along the device channel. Then the 
diagonal components of the resistivity tensors arc con- 
tinuous across the DW, and the off-diagonal components 
change sign. For small /3 the longitudinal current com- 
ponent Jx will be dominated by the uniform background 
current j, except within a distance ~ Pw of the sides of 
the device near the DW where the current perturbation 
is concentrated. Neglecting this edge correction, it then 



3 




FIG. 2: Calculated AMR contribution to the resistance across 
a DW for different values of the anisotropy /3. Symbols, nu- 
merically calculated values; solid line, values using Eqn. 
A film thickness t = 100 nm, resistivity p = 3 x 10^* Qm, 
misalignment 5 = 0°, DW angle i9 = 0° and magnetisation 
angle (p = 45° have been used. Insets: the effect of varying 
the DW angle i9 and misalignment 5. 

follows from Eqn. ^ that immediately on either side of 
the DW 

Jj,(a;o±0+,y) ^±^sin2(^. (10) 

Ignoring anisotropy, the slowest decaying current pertur- 
bations decay like^^ exp— 7r|a;|/ii; and by assuming that 
Jy decays like this from the interface value we can 
evaluate the integral in ^ to get for the AMR contribu- 
tion to the average longitudinal resistance 

Ramr = -Rdtt- sin^ 2(/3 (11) 

ZTT 

where i?n = p/t is the sheet resistance. This result shows 
how the circulating currents give rise to a negative con- 
tribution to the resistance across the DW. 

We have also performed numerical studies of the cur- 
rent distributions and resulting fields and voltages in the 
presence of DWs. The numerical solution is not restricted 
to the specific configuration that was assumed in deriving 
the analytic estimate for i?AMR, so as well as enabling an 
assessment of the accuracy of this expression obtained 
we are also able to include the effects of misalignment 
of the hard axis with respect to the device channel, and 
the angle of the domain wall. Some results are given in 
Fig. O The solution is obtained by introducing a stream 
fimction ip{x^y) that is related to the current density 
via J = (d4'/dy,—d'ijj/dx), thereby ensuring that cur- 
rent continuity Eqn. (j2a[) is satisfied. Combining Eqns 
(|2b[) and (Pc|) then results in a non-separable elliptic par- 
tial differential equation^^ that we solve for ip via the 
multigrid relaxation method. 

In Fig. [3^ we show typical results from our numerical 
studies, displaying the variation in the longitudinal resis- 
tance as a DW inclined at an angle = 20° passes along 
a device channel of width w = 30/im and with voltage 
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FIG. 3: (a) Variation in average resistance along the sides 
of the device channel (between a; = and x — 180 fim) as 
a function of the position of the DW. The aspect ratio is 
l/w — 6 and the DW inclined at i9 = 20°. See the text for the 
other parameters used. The inset shows the induced current 
flow, (b) Similar to (a) but for different material parameters 
(see text) and dimensions: the resistance is measured between 
I = and X = 200 fim, l/w = 2 and i9 = 50°. 

probes separated by I — 180/im. The average resistance 
of the two uniform magnetisation states has been sub- 
tracted: ARxx = (R) — {elx + Qxx)^ / C^^t) . A general 
linear variation in the resistance is seen, except when the 
probes are within the range of the circulation currents 
induced by the DW; these cause a rapid variation over 
a distance w tan i9 -I- 2w /tt as expected from geometri- 
cal considerations and the discussion following Eq. (jlOp . 
Furthermore, we see that the calculated resistance lies be- 
low a straight line interpolation performed between the 
two asymptotic channel resistances of the two magneti- 
sation states. In this calculation we use values that cor- 
respond as best as possible to the system reported in Fig. 
4 of Ref. [ll|: film resistivity g = 3x 10~* Jim, thickness 
t = 100 nm, and anisotropy-^ f3 = —0.03. The magneti- 
sation orientations within the two domains are taken to 
he — S + (j), (j)2 = S — (j) where 4> = 37° due to uniax- 
ial anisotropyi^iii and the misalignment S = —0.28° (the 
difference between the asymptotic resistances, 

(eL - = ■^sin2(?!)sin25, (12) 

wt wt 

is then 5 $1 as found in Ref. [HI). Using these values, 
numerically we find the resistance is lowered by 0.33 f2 
as a result of the eddy currents induced by the DW. 

In Fig. [2] we compare numerical values for i?AMR found 
in a number of similar calculations to that just described, 
with those obtained using Eqn. ()11|) . The numerical re- 
sults also display the 0^ dependence and, as expected 
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given the approximations made in estimating the current 
integral, our analytic expression overestimates the actual 
resistance, we find by approximately 15% when -d — 0. 
This value of i?AMR is further reduced at DWs inclined 
relative to the current direction, but is relatively insen- 
sitive to the misalignment angle (insets in Fig. [2|). Thus 
Eqn. pTjl has some value in estimating the AMR contri- 
bution to the DW resistance, but numerical calculations 
are required for accurate estimates. 

In Fig. we show the calculated longitudinal resis- 
tance for a second case, with parameters chosen to corre- 
spond to the device reported in Fig. 1 in Ref. TT]. This 
was the initial device studied experimentally, in which 
the misalignment of the hard axis with the device chan- 
nel is greater. We use g = 4 x 10""* ilm, t = 150 nm, 
P = —0.03, with the channel width w = 100 /im and 
voltage probes separated by 200 /im. Also, — —37°, 
6 — 1.5°. The greater structure exhibited by AR^x in 
this case is due to a larger DW inclination {d = 50°) and 
the smaller aspect ratio Ijw oi the device, and results in 
a linear variation over only a short range of DW positions 
midway between the voltage probes. The precise results 
are rather sensitive to the value of "d. However, generally 
we find that in the linear region the large spatial extent of 
the eddy currents still affects the slope of the resistance 
curve, which no longer coincides with a linear interpo- 
lation of the asymptotic resistances between the probe 
positions. The dashed line in Fig. which is parallel to 
the linear section of the resistance curve connects points 
some 30% further apart than the probes. The calculated 
resistance is again lowered due to the AMR, but the dif- 
ference of —0.16 ri is smaller than the value (—0.18 i7) 
found if the distance between the voltage probes is in- 
creased so as to fully contain the eddy currents. 

The magnetisation profile within the wall can also 
give rise to a negative AMR. However, the contribu- 



tion we describe above dominates here. Assuming a 90° 
Neel like wall with magnetisation rotating like (f){x) — 
— (1/2) tan^^ sinhx/A, where ttA is the wall width, gives 
a contribution to leading order of A/(— 2/3w)i?AMR, or 
just a few percent of the contribution from the circula- 
tion currents. Other wall profiles in which the spin ro- 
tates out of the plane lead to the same conclusion. Only 
if the DWs in this system were 180° walls would the in- 
wall contribution be significant, since then the circulation 
current contribution (jlip vanishes. 

Comparing with experiment, the DW resistances re- 
ported in ll| for the devices modelled in Figs [5^ and 
are — 1.0±0.2 and — 0.44±0.5 Vt respectively; a third set 
of devices similar to that of Fig. [3^ but with w = 60 /im 
gave —0.3 ± 0.2 Vl. The corresponding i?AMR values we 
find are —0.33 il, —0.16 and —0.33 fi. The previously 
neglected AMR contributions to the resistance across the 
DW make a major contribution to, and can largely ex- 
plain, the negative values observed, with the exception 
of one set of devices where a true negative intrinsic DW 
resistance may indeed have been observed. Clearly fur- 
ther experiments are required to clarify the situation, 
before attempts to quantitatively account for the DW 
resistance's, can be properly assessed. For these, devices 
with a large aspect ratio l/w, and containing DWs orien- 
tated normal to the device channel, are clearly desirable. 

To summarise, we have identified a significant 
anisotropic magnetoresistance contribution to the neg- 
ative domain wall resistivities recently observed in mi- 
crodevices fabricated from (Ga,Mn)As epilayers. We de- 
rive an analytic estimate of the magtitude of this contri- 
bution, and report calculations of the channel resistance 
as a DW is moved through the device which provide a 
good description of the experiments. 
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